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Abstract 

We introduce the A-D-E resonance factorized models as an appropriate analytical 
continuation of the Toda S-matrices to the complex values of their coupling con- 
stant. An investigation of the associated Casimir energy, via the thermodynamic 
Bethe ansatz, reveals a rich pattern of renormalization group trajectories interpolat- 
ing between the central charges of the Gi Gk/Gk+i GKO coset models. We have 
also constructed the simplest resonance factorized model satisfying the "i?i>'^" -property. 
From this resonance scattering, we predict new flows in non-unitary minimal models. 
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1 Introduction 



Conformal invariance has played an important role in our understanding of the properties 
of two-dimensional systems at the critical point |l|, Away from the critical point, the 
corresponding field theory can be considered as a conformal model perturbed by an ap- 
propriated relevant operator. The important insight added by A. Zamolodchikov was 
that certain perturbations may possess infinitely many non-trivial conserved charges, hence 
integrable even in the off-critical regime. This observation allowed him to construct the 
on-mass-shell solution, i.e. the mass spectrum and the exact S-matrices, of the Ising model 
in a magnetic field 0. In this context, Al. Zamolodchikov has pointed out that by 
applying the thermodynamic Bethe ansatz (TBA) one is able to recover the ultraviolet 
properties of the ground state from the respective S-matrices. In this sense, the TBA makes 
an interesting bridge between the properties of the infrared and the ultraviolet regimes. In 
principle, this technique can be applied to any factorized scattering theory, even for those 
that a priori lack a background conformal field theory. An a posteriori analysis of the be- 
haviour of the finite-size corrections to the ground state may reveal amazing features not 
previously known. Recently, one example has been put forward by Al. Zamolodchikov 
in the case of a resonance S-matrix consisting of a single particle scattering through the 
following amphtude, 

era a \ - sinh(g) -2cosh(go) 
'^^^'^°^-sinh(^)+^cosh(^o) 
where is the resonance parameter. 

Let us briefly summarize Al. Zamolodchikov's results 0. Although the ultraviolet limit 

associated with this S-matrix is always governed by a theory with central charge c = 1, 

at intermediate distances (for Oq » 1) the pattern of the renormalization group (RG) 

trajectories is surprisingly rich. The RG trajectories evolve to the infrared regime by first 

interpolating between the central charges Cp = 1 — 6/p{p+ 1) of the minimal models. Each 

fixed point Cp takes approximately y of the RG "time" before crossing over to the next 

lower critical point Cp_i. At this point we recall that such crossover phenomena can be 
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induced by the perturbation of the least relevant field operator 0i 3 of the minimal models 



[|, |T0[- Recently, it has also been argued in ref. [|Tl|] that the field 0i,3 can produce 
crossover behaviour even in non-unitary minimal models. 

The purpose of this paper [] is to present and analyze a large class of resonance factorized 
scattering theories based on the simply laced Lie algebra G. These S-matrices will define 
a family of one parameter field theories G{6o), depending on the resonance parameter ^o- 
We point out a connection between these resonance S-matrices and those from the A-D-E 



Toda field theories [|T6|, |T^, |T^, through an appropriate analytical continuation of 
the Toda coupling constant to complex values. This allows us to express the resonance 
S-matrices amplitudes in a rather simple and useful form. We study the Casimir energy 
Eq{R, 6q) at moderate distances R via the TBA approach. The behaviour of the associated 
RG trajectories are then related to fiows of certain deformed GKO coset models |]T2 



Finally, we propose a simplest resonance factorized model satisfying the "0^" -property. We 
argue that the respective RG trajectories can be related to new fiows in the non-unitary 
minimal models. 

The paper is organized as follows. In sect. 2 we introduce the A-D-E resonance models 
and write down the respective TBA equations. Sect. 3 is devoted to our numerical checks of 
the behaviour of the ground state energy Eq{R, 6q) for finite values of the volume R. Sect. 4 
describes the "0^" -resonance model and our arguments to predict new fiows in non-unitary 
minimal models. Sect. 5 contains a discussion of our results. In appendix A we summarize 
some useful relations in order to compute the logarithmic corrections to the Casimir energy 
when i? — > 0. 



brief account of this paper has recently been presented in ref. | 
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2 The A-D-E resonance S-matrices and the TBA equa- 
tions 



We start this section by describing the basic properties of the resonance scattering that we 
have constructed. We assume that the S-matrices are diagonal and factorizable, i.e., there 
is no particle production. In this case, the only constraints are due to the unitarity and 
crossing conditions of the two-body amplitude Sa,b, Q 

SaMSa,b{-0) = l, Sa,biO) = SbA^n - 9) (2) 

where a stands for an antiparticle. 

We also assume that the resonance amplitude Sa,b{(^, Oq) can be written as a product of 
a background amplitude S'^l"^{9) and the resonance part Sl^^{9,9o), namely 5'a_b(6', 6'o) = 
S^b"'i^)Sa'^bi.^y ^o)- The physical meaning of this hypothesis is that the background and the 
resonance interactions appear in different spatial regions, and therefore the phase-shifts are 
simply added. More complicated cases will involve correlations that will not be considered 
in this paper. However, as we have already stated, even this simple ansatz will lead us to a 
rather rich pattern of RG trajectories appearing in deformed conformal field theory. 

In order to introduce the A-D-E structure we assume that the particle spectrum and the 
background S-matrices, S^l^{9), are given by the so-called minimal solution of Eq.(2). For 
details and respective properties see ref. |]16|, [l^, |19|, |2^. The simplest solution of Eq.(2) 



to the resonance amplitude Sl^^{9, 9q), invariant by the symmetry ~^ ~^0) appears in the 
case of the Ajv-i Lie algebra. The S-matrix for the fundamental particle is given by 

oresfn n ^ ^ siuh ^ (g - - siuh |(g + - ) 

^ ' °^ sinh \{9-9^ + ) sinh \{9 + 9^ + ij^) ' ^""^ 
where the resonance poles are located at = ±0p —iff- The other amplitudes can be 
computed by a straightforward bootstrap. For N=2, we recover Al. Zamolodchikov's S- 
matrix, Eq.(l). 

^In our case the antiparticle appears as a pole in the particle-particle scattering amplitude, and therefore 
the reflection amplitude is null in order to satisfy factorizability [15|. 
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Similar expressions can be found for the D and E groups, by consistently solving Eq.(2) 
in analogous manner as pursued in refs. [|16|, 0, |T8|, |T9|, However, we prefer to take 
a short route. First, we notice that the amplitude Sl^i{6,6o) above is related in a simple 
way to the so-called Z-factors of the Ai\i_i Toda field theory ||16|. These factors encode the 



dependence of the coupling constant a present in the Lagrangean, through a function b{a) 
[0. Then, by simply taking b{a) = ±i9o + in the Toda Zi 1(6', 6(Q;))-factor we recover 
the amplitude of Eq.(3) 0. In the general case, we have verified that the resonance S-matrix 
satisfying the properties mentioned above is given by 0, 

sii^ie, e,) = Za,,{e, h{a) = ±ie, + ^) (4) 



where Zafi{0,b(a)) are the Z-factors appearing in the A-D-E Toda field theory |]TB|, [T7 



p!9| , |2(| and h is the Coxeter number. 

Using Eq.(4) we can write down a closed expression for all the Sa,biO,Oo) amplitudes. 
First it is convenient to define the functions, 

i^aAO,Oo) = ~^-^\ogSa,tie,(^o) (5) 

which will also be important in the TEA equations. In terms of its Fourier component 
^a,fe(fc,6'o) defined by, 

_ roo 

i^a,b{k,Oo)= / e^'Va,6(^^,^o) (6) 

J —oo 

we find the remarkable matrix identity. 



^,.,cosh[f]-/,V2 

cosh[^] - cos{k9o) ^ ' 



where lafi = 2 — Ca^ is the incident matrix of the G=A,D,E Lie algebra. Finally, the 
amplitude Sa,b{(^, Oq) can be written in terms of a standard integral representation, 

SaAo.eo) = exp J^^4,^^,{k,eof-^^^dk] (8) 



■^The case N=2 was first noticed by Al. Zamolodchikov in the context with the sinh-Gordon modeL 
^Here we have followed the notation of ref. p6| . There are other similar notations in the literature 

^, |l^, For example, in the case of ref. ,by taking B{a) = 1 ± iOq^ we achieve the same results. 



We now start to discuss the TBA equations and their properties. The TBA approach 
allows one to compute the Casimir energy Eo{R, 9q) on a radius of length R. The ground 
state energy is written in terms of pseudoenergies ea{0), and for diagonal S-matrices there 
is one such function ea{0) for each stable mass of the spectrum. The Eq{R,6) is given 
by, 



N 



27r 



a=l 



Y^nia ddcosh{d)La{e) 



(9) 



where La{9) = ln(l + e"'^"^^-'), and N is the number of particles. 

The volume dependence is encoded by a set of integral equations for the pseudoenergies 
ea{0), depending on the kernel ipa,b{0,(^o) as follows, 

N 



1 roo 



(10) 



where z/a(6') = ma/2cosh(^) and in our case rria are the mass gaps of the A-D-E spectrum 
[0, |T^, |T^, It is possible to rewrite Eq.(lO) in a more suggestive way, adopting an 
approach similar to that of ref. |21|. Using Eq.(7) and the following identity [pi], 



N 

we then obtain a set of functional equations for the functions Ya{6) = e'^°-^^\ 

-1 



(11) 



beG 



1 + K 



-1, 



-1 



(12) 



Motived by Al. Zamolodchikov's discussion of the Ai it is possible to rewrite Eq.(12) 
as an interesting discrete equation if we restrict the rapidities 6 on the lattice variables 



m,n 



9 + mOQ + ^n. Using the identification yj 



= Ya[era,n) WC haVC, 



-ym,n+l-\^iri,n—l 
a a 



m,nila.b 



n 



1 



1 a 



(13) 



where b is an index characterizing the nodes of the Dynkin diagram of G=A,D,E Lie algebra 
and c is the same index for the A Lie algebra. 

Eq.(13) is defined on the product space of the Lie algebras G (A,D,E) and A, namely 
G ® A. Eq.(13) also assumes a form similar to those obtained from certain perturbed GKO 



coset models p4|, It is worth mentioning that a similar equation can also be defined 
in the dual space A ^ G, provide that we take Ya — > Y^^ in Eq.(12). Therefore they are 
self-dual when G = A. In this sense these two pairs of equations are richer than the TBA 
equation that we start with, Eq.(lO). However, the TBA equations can be easily recovered 
from Eq.(12) by specifying the initial conditions of the functions Vaid)- 

In the next section, we show that even the simple guess for the energies nodes Ua, i.e. 
Va{0) = maRcoshlO), will generate a rather interesting pattern of RG trajectories. It turns 
out that these trajectories resemble those of unstable particles decaying to intermediate 
states until reaching the vacuum. At each intermediate step the particle has a life-time of 
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(the RG "time") before decaying again to the next state. 



3 The numerical work 

We begin by analyzing the ultraviolet behaviour of Eqs.(9,10). Standard manipulations of 
the TBA equations p give us the following behaviour of the ground state energy, 

vrr 

Eo{R,9o)c^-—,R--0 (14) 

where r is the rank of the respective Lie algebra. 

From Eq.(14) we conclude that the background conformal field theory has central charge 
c = r |22[. In appendix A we discuss the next to the leading order in R corrections to 
Eo[R,6o). The first correction is logarithmic in R , in contrast with the fractional power 



behaviour of the typical non-resonance factorized theories |^T|. In terms of the function 
c{R,eo) = -^Eo{R,eo) our result reads 

c{R,eo)=r-^-^^^±i^j:C-i (15) 

a,b 

where X = In(^). 

The behaviour of the function c{R, 6q) becomes particularly interesting at intermediate 
lengths of R. In this case one needs to solve numerically the integral equations (10) using, 
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e.g., standard iterative procedures. We start by analyzing the simplest case of our proposed 
resonance S-matrices, namely the A2 theory ^ In fig. l(a,b,c,d) we show the behaviour of 
the function c{R, 9q) for = 0, 10, 20, 40. c(-R, 0) behaves as a smooth function interpolating 
from the ultraviolet regime, c{R — * 0, 0) = 2, to the infrared region. On the contrary, by 
increasing 6'o, it is clear that c(-R, 6*0) starts to form plateaux around the following fixed 
points, Cp = 2(1 — 12/p{p+ l)),p = 4,5,6,... . These critical points correspond to the 



conformal field theories possessing a Z(3)-symmetry |2^. Since the numbers of plateaux 
increase with 6*0 one expects that when 6*0 — > 00, all the fixed points will be visited by the 
RG trajectories until finally reaching the infrared region. 

It is also fruitful to examine the same pattern in terms of the /3-function along the 
RG trajectory. In this case, the plateaux will correspond to zeros of the beta function. 
Borrowing Al.Zamolodchikov's definition |^ of the /3-function, we have 

m = -^ciR,eo), g = r-ciR,9o) (16) 

where the "coupling constant" g has been normalized to zero in the ultraviolet regime. 

In Fig. 2(a,b,c) we show P{g) for 9q = 10, 20, 40. One can notice that the zeros of P{g) are 
precisely at the points g = 24/p(p+l),p = 4, 5, 6... By increasing 6*0, more zeros are going to 
appear in l3{g), however keeping the general shape of l3{g) for those that have already been 
formed. We stress that our findings are very much in agreement with Al. Zamolodchikov's 
results for the Ai case. In general, for the proposed A-D-E scattering theories, the function 
c{R, 6q) will form plateaux around the values of the central charge of the Gi^Gp^h/Gp-h+i 
GKO |]T^ coset construction, namely Cp = r(l — h{h + l)/p{p + l)),p = h + 1, h + 2, ... . 
In Fig. 3(a) we illustrate this pattern for the first plateaux in the models A2, A^, D4, which 
share some common central charges. In Fig. 3(b) we plot the respective f3{g) function. 

What we can learn from this numerical computation is as follows. Each time that 
X ~ —{p — h)^ the function c{R,9q) crosses over from its value to the next (up) value 
dpj^i Q- This indicates that at these values of X the TBA equations may present a rather 
^ A similar study has been performed by Al. Zamolodchikov in the case of ^1 

^ From our figures 1(c), 1(d), 3(a) we notice that the "RG time" ^ accounts for the plateau and for 



special behaviour. Indeed, the hnearization of these equations around X ~ — (p — /i)^ (or 
better the functional Eq.(12)) gives us a set of 6'o-independent equations [], describing the 
RG flows in the Gi ® Gp-h/Gp-h+i coset models perturbed (positive perturbation) by an 
operator $p with conformal weight A^^ = 1 — h/{p + 1) |]25| , ^]. On the other hand, 
the infrared finite size corrections will be dominated by the dual field $p with conformal 
dimension A^^ = l + /i/(p — 1)0. Therefore, for each plateau "p" the ultraviolet behaviour is 
governed by $p and finally attracted to the infrared region by $p. It is tempting to propose 
that such behaviour may be reproduced if one considers the critical theory perturbed by 
the combination X^p + X^p (A > 0) P, In fact, in the case of the minimal models. 



Lassig has identified such a phase for A < 0, by using perturbative RG calculations 

[|, p. The generalization of this approach to all A-D-E resonance models are technically 

complicated, because the necessary set of structure constants of the cosets models are not (up 

to now) explicitly known. However, we point out that the field $p has the correct conformal 

dimensional to satisfy some of the arguments used in ref. ||2^. Moreover, at least for the 

ground state, our numerical results are consistent with this picture. 

Finally, we remark on some numerical observations concerning the behaviour of the 

functions La{6). In fig. 4 we show this behaviour for the model. The Z(4)-symmetry 

assures that Li{6) = L^{9) and therefore we only plot the functions Li{9) and L2{9) for 

> {La{6) are even functions). Roughly it seems that these functions differ only by an 

overall constant (for each "p" ) . This fact may be expected considering that these functions 

are smooth and should assume constant values at the plateaux. In fig. 5 we have plotted 

the first derivative (in the variable 6) of Li{6) and L2{6) for X = —70. We believe that this 

the finite size corrections, and both are model dependent. For example, in the case of the Ajv-i a rough 

estimate of the plateau's length is 

^ One should notice that, e.g., the kernel V'l.il^i^o) has deep peaks a.t 9 ~ ±9q,0, appearing as the 

interacting (A Lie algebra structure) and the diagonal (G Lie algebra structure) of the linearized TBA 

equations. 

* For p = h + 1 this field is not present in the Kac table and it is replaced by the spinless combination 
TT of the stress energy tensor T . 
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result may motivate us to derive more exact results from the TBA equation (10). Anyhow, 
for the moment, this information has been used (Appendix A) to estimate the logarithmic 
correction to the Casimir energy in the ultraviolet limit. 



4 The "(/)^" resonance S-matrix 



Let us consider a solution of Eq. (2) possessing a single particle a. We also assume that the 
pole structure in the scattering amplitude produces the particle a itself. Besides Eq.(2) the 
S-matrix Sa,a should satisfy the "0^" bootstrap condition, 

SaAd^ do) = SaA0 + Oo)SaAO ' ^o) (17) 

Following the considerations of section 2 , the simplest solution of Eqs. (2,17) are given 

by, 

tanh |(g + zf ) tanh |(g - - ) tanh |(g + - zf ) 
"'"^ ' °' tanh|(^-zf ) tanh 1(0- 0p + if) tanh i(^ + ^p + zf) ^ ' 

The background part of SaA^j^o) is the S-matrix of the perturbed Yang-Lee edge 
singularity pSf. The Toda related field theory is that one analyzed by Arinshtein et al 



[1^ and known as the Shabat-Mikhailov model. This theory has only one field and the 
interaction is through the potential V{(f)) = e°"^ + e"^""^, a being the coupling constant. 
As has been shown by Mikhailov p9|, this model is integrable and can be considered as a 



"reduction" of the A2 Toda field theory provided that the two fields (pi, 02 in the A2 model 
satisfy the relation (pi = —(p2 = (p (see ref. |]2D|). From the S-matrix point of view, this can 
be seen by defining the amplitude SaA^^ ^0) as 0, 

SaAG^o,) = SiAo,e,)SiAe,eo) (i9) 

where S'i^i(6', 6'p) and 5*1^2(^,^0) are the A2 amplitudes. 

It follows from Eq.(19) and Eq.(lO) that the pseudoenergie (e(6')) of the model defined by 

Eq.(18) and those from the A2 Toda theory (ei(6'), €2(6*)) are equal, namely e{6) = ei{9) = 
^ A similar reduction as in Eq.(19) is possible for any A^v-i theory, N odd. The next task is, however, 
the identification of the reduced model for any kind of crossover behaviour in non-unitary models. 
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e2{0)- Hence the ground state energy associated to the S-matrix of Eq.(18) is precisely half 
of that of the A2 model analyzed in sect. 3 . Therefore, the plateaux will now form around 
the values Cp = 1 — 12 /p{p + = 4, 5, .. . At this point it is important to recall that the 
Casimir energy in non-unitary minimal models Mr depends on both the central charge and 
the lowest conformal dimension A^m ( effective central charge) pO|, 



Eo --^(-24A„i„ + c) ~ --^Cef (20) 
QR 6R 

where Cg/ = 1 — 6/{pq)- 

This fact strongly suggests that we are dealing with non-unitary minimal models. Indeed, 
comparing Cg/ with Cp = 1 — 12/p{p +l),p = 4, 5, ... we identify the following series of non- 
unitary minimal models: (q=p/2, p even) and M g+i (q=(p-l)/2, p odd). We also 
recall that for these models the field 0i^2 has the minimal conformal dimension. Motived 
by these analogies with the A2 resonance factorized theory it is also interesting to find the 
fields that may drive the system to the staircase pattern of section 3. From the expected 
finite size corrections to the ground state we identify the fields 02,i and 0i,5 as those that can 
induce these models to the crossover phenomenon. More precisely, we expect the following 



behaviour for M g+i , 

2q+l 



and for M q , 

2q+l 



Mq+i +(/)9i -> g = 2,3,... (21) 

2q + l ' 2q + l 



M^+(/)i.5^M g = 3,4,... (22) 

29+1 2{9-l) + l 

Moreover, the infrared regime is also governed by the operators (f)2^i and (pi^^. An excep- 
tion is the flow M3/7 + 01 5 —>■ M3/5 , where the field 0i^5 is not present in the M3/5 Kac-table. 
In this case we think that 0i,5 should be replaced by the level 2 descendent of 0i,3. To the 
best of our knowledge this pattern in which the fields 02,i and 0i^5 interchange their roles as 
relevant and irrelevant operators along the RG trajectory is new in the literature. In fig. (6) 
we show this predicted behaviour R. 



It is notable that the central charge increases only in the flow defined by Eq.(22). Since these models 
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In order to give more support to this prediction we have numerically studied the spectrum 
of the Hamiltonian of the simplest flow, namely M3/5 + 02, 1 ^^2/5- We have constructed 
the Hamiltonian of this perturbed conformal field theory on a torus of radius R using the 



truncated conformal approach ||3^, |3^, In fig. (7) we show the evolution of the first 
40 eigenvalues as a function of the volume R. The 02,i perturbation allows us to divide the 
Hilbert space in two sectors of the fields [/, 02, 1] (dashed lines) [0i,2, 0i,3] (solid lines) and 
their descendents. The striking feature is the absence of any crossing levels, that makes it 
difficult to interpret fig. (7) as a consistent integrable massive theory. We recall that in an 
integrable massive theory, there will be an abundance of level crossings between states that 



are distinguished by their particle content |33, Bq]. For example, this will be the case for 



momentum lines accumulating in the lowest particle threshold and other n-particle states 
present in the theory. Moreover an a posteriori analysis of the (possible) mass gaps are not 
consistent (within our numerical precision) with, e.g., the exponential split of the ground 
state. Although, we are not able to estimate the infrared region using this approach, fig. (7) 
has a rather remarkable resemblance to typical massless crossover spectrum. We believe 
that these observations bring some extra support to the prediction of Eqs. (21,22). Finally, 
some additional comments are in order. 

A first remark concerns a possible Ginsburg-Landau |^ interpretation of the simplest 
fiow M3/5 + 02,1 M2/5. The fusion algebra for the operator 0i,2 with the lowest dimension 
in the M3/5 model is. 



'1,2 ^ (Pl,2 ~ J (Pl,3, 91,2 ^ (Pl,3 ~ 92,1 



+ 01,2; 01,2 ® 02,1 ~ 01,3 (23) 



By identifying the operator 0i 2 with the elementary field $ of the Ginsburg-Landau 
theory, it follows from Eq.(23) that 0i,3 = and 02,i = The first descendent L_iL_i0i_2 
being an odd field should be identified with The Ginsburg-Landau potential consistent 
with this last identification is V{^) = which is in accordance with the Z{2) symmetry of 
are non-unitary they do not present any contradiction with Zamolodchikov's c-theorem For other 

interesting behaviour in polymer systems see the discussion in ref. p7| . 
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M3/5 model. Therefore, perturbing with the field 02, 1 one can drive the M3/5 model to a $^ 
class of universality, which is precisely that of the Yang-Lee edge singularity [^] It seems 
important to find the Ginsburg-Landau picture of other non-unitary models. However, we 
have noticed that in models with several negative and positive dimensions it seems quite 
difficult to find a consistent description in terms of only one basic field $. 

A second observation concerns the ambiguity of sign of the perturbing coupling constant 
in the case of the 0i,5 operator. The convention that is assumed in unitary models is that the 
positive sign is related to the massless crossover and the negative sign to a massive behaviour. 
Since these non-unitary models are somewhat related (ground state) to an unitary theory 
{A2 model) one may expect that the same picture holds. However, in the case of non-unitary 
theories we may also have phases in which the system starts to present complex eigenvalues 
after reaching a certain value of R . We believe that a further investigation is needed to 
resolve this issue, which may end up in discovering new S-matrices for non-unitary models. 



5 Conclusions 

We have studied resonance factorized models based on the A-D-E Lie algebra. By using the 
thermodynamic Bethe ansatz we have analyzed the respective Casimir energy on a radius 
of length R. At moderate distances of R we find a rich pattern of RG trajectories associated 
with typical fiows in the deformed GKO coset models. We have pointed out a connection 
between the A-D-E resonance scattering and one from the Toda field theory when its cou- 
pling constant assumes special complex values. However, the meaning of this analytical 
continuation in terms of a well defined background quantum field theory is still to be found. 
Moreover, it would be interesting to investigate whether this is the only possible way leading 
to factorized resonance S-matrices which can be interpreted as a deformed conformal field 

theory. We have obtained a general functional equation for the functions Ya{0) = e*^^^-* and 
We also recall that the identification of the operator 0i 3 is in reasonable agreement with a Z{2) 
spontaneously broken symmetry (Ising-like) found in its spectrum p^. 
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one may start with this equation instead of considering the A-D-E resonance S-matrices. 
We recall that similar functional equations appear in integrable lattice models in the context 
of the inversion relations to the transfer-matrix . It seems also that these equations hide 
a large amount of information yet to be explored ^ . 



The construction of a resonance factorized S-matrix satisfying the "0^" -property has led 
us to predict new flows in the non-unitary minimal models. In the case of the simplest 
crossover, the spectrum presented in section 4 may give us a guess on how the fields in 
the ultraviolet regime will evolve to the infrared region. However, in the more general 
cases considered in this paper, this property has still to be disentangled. It would be also 
important to identify the lattice models possessing such properties, the yl2-RS0S models 
being a probable guess [HO . 
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6 Appendix A 



This appendix is concerned with the logarithmic corrections appearing in Eq.(15). Our 
analysis is fairly parallel with that of ref.0, in order to generalize it to the case of many 
pseudoenergies. First we notice that Eq.(7) permits us to write the expansion 

where up to order k"^ we have, 



^S=^a,.; V;S = 2C-^J^o' + ^) (25) 



Using this expansion in Eq.(7), in the regime X — oo, we find 



OO „/",(2") 



m.e^ + ln(l - e'^^^^)) = E E (26) 

b n=l l^'^J- 

where Lf"^) = f^eLb{0). 

Following ref. [0, the central charge in the region X ^ y, 6 — 2X » 1, y << 0, satisfies 
the differential equation, 

^^c{X,y) = -J2m,eyLa{y) (27) 
Q dy V 

and one possible ansatz for its solution is given by, 



oo ^/,(2n) 2n-l 



" ^ a,b n=l l^^J- k=l a -^0 a 

(28) 

The leading logarithmic correction is obtained by expanding ln(l — e~*) = — e^* + 
and by ignoring the term in Eq.(26). The final result is, 

c{X, y) = r + ^ ^2^«i.?^ - ^ E (29) 



provided that the La{0) satisfy, 

E^4b^^(^) + e-'"^'^ = o (30) 



^ 2 
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Based on our numerical results, it is a reasonable approximation to consider that the 
first derivatives of the La{9) are index a independent. Therefore, we can have the following 
ansatz as a solution of Eq.(30), 



Lai9) = In 



sin^ - 7) 



(31) 



where A is fixed by periodicity to be A = 2(x-7o) ' ^ constant (not determined by this 

approach) such that when X ^ 00, j ^ Jq. Substituting Eq.(31) in Eq(29) we have, 

■J /^a,b 'ra,b -t^ /oo\ 

c - ^ - 2 ' X^oo (32) 
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Figure Captions 



Fig. l(a,b,c,d) The scaling function c{R, Oq) for four values of the resonance parameter 
^0- (a)^o = 0, (b) 00 = 10, (c) 00 = 20 and (d) Oq = 40. 

Fig. 2(a,b,c) The beta function P{g) for three values of the resonance parameter ag- (a) 
9o = 10, (b) 9o = 20 and (c) 9o = 40. 

Fig. 3(a) The first RG plateaux for the models A2 (short dashed hue), (solid hne) and 
D4 (long dashed hne). 

Fig. 3(b) The beta-function (3{g) for the models A2 (short dashed line), A3 (solid line) and 
D4 (long dashed hne). 

Fig.4 The functions Li{9) (solid hne) and L2{9) (dashed line) of the A3 model for ^0 = 40 
and X = -70. 

Fig. 5 The first derivative of the functions Li{9) (solid line) and L2{9) (dashed line) of the 
^3 model for ^0 = 40 and X = -70. 

Fig. 6 The flow pattern in the non-unitary minimal models M 1 ; M g+i , g = 2, 3, .... The 

2<3+l 213+1 

horizontal(vertical) arrows represent the relevant (irrelevant) operators defining the ultravi- 
olet (infrared) corrections to the fixed point. 

Fig. 7 The first 40 levels of the spectrum of the model M3/5 -|- 02,i- The dashed (solid) hues 
correspond to the sector [/, 02,i] ([0i,2, 0i,3]) and its descendents. 



18 



